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Abstract 

o 

Burke's theorem can be seen as a fixed-point result for an exponential single-server 
$H queue; when the arrival process is Poisson, the departure process has the same distribu- 

^2 tion as the arrival process. We consider extensions of this result to multi-type queues, in 

<^ which different types of customer have different levels of priority. We work with a model 

l/^ of a queueing server which includes discrete-time and continuous-time M/M/1 queues 

as well as queues with exponential or geometric service batches occurring in discrete 
time or at points of a Poisson process. The fixed-point results are proved using inter- 
Ql^ changeability properties for queues in tandem, which have previously been established for 

Qh one- type M/M/1 systems. Some of the fixed-point results have previously been derived 

as a consequence of the construction of stationary distributions for multi-type interact- 
ing particle systems, and we explain the links between the two frameworks. The fixed 
points have interesting "clustering" properties for lower-priority customers. An extreme 
I I case is an example of a Brownian queue, in which lower-priority work only occurs at a 

, set of times of measure (and corresponds to a local time process for the queue-length 

^ process of higher priority work). 

1 Introduction 

cn 

One of the most famous results in queueing theory is Burke 's theorem [7 . Consider a queue 
in which available services occur as a Poisson process of rate /i (a so-called ./M(/i) /I queueing 
server). If the arrival process is a Poisson process of rate X < ^ (independent of the service 
I> process), then the departure process is also a Poisson process of rate A. We may say that the 

arrival process is a fixed point for the server. 

In this paper we consider the question of fixed points for queues with two or more classes 
of customer (with different levels of priority). When a service occurs in such a queue, it is 
used by a customer whose priority is highest out of those currently present in the queue. We 
will see that a two-type fixed point can be constructed using the output processes (consisting 
of departures and unused services) from a one-type queue. Then in a recursive way, a fixed 
point with m > 3 classes of customer can be constructed using the output of a queue whose 
arrival process is itself a fixed point with m — 1 classes. 

Except in the familiar one-type case, the fixed points are not Markovian. In particular, 
one observes clustering of the lower-priority customers. 

In the paper we work with a queueing model which is somewhat more general than the 
./M/1 queue described above. Our basic model is of a discrete-time queue with batch arrivals 
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and services. Let 5„ be the amount of service offered at time n. We obtain fixed-point results 
for the case where Sn are i.i.d. and each 5„ has so-cahed "BernouUi-geometric" distribution, 
i.e. is equal to the product of a geometric random variable and an independent Bernoulli 
random variable. By taking appropriate limits where necessary, this model covers a variety of 
previously considered queueing servers, for example discrete-time -/M/l queues [17], discrete- 
time queues with geometric or exponential service batches [SI [HI HZ], continuous-time ./M/l 
queues as described above, continuous-time queues with geometric or exponential service 
batches occurring at times of a Poisson process, and Brownian queues [13 [2S]. Versions 
of Burke's theorem and related reversibility results were proved for this Bernoulli-geometric 
model in [21) . 

Some such fixed-point processes were already constructed in certain cases (M/M/l queues 
in continuous or discrete time) in [12] in the context of stationary distributions for certain 
multi-type interacting particle systems. In this paper we give a more direct proof of the 
fixed-point property, which relies on properties of interchangeability for queueing servers. 
Weber [3D] showed that for a tandem queueing system consisting of two independent ./M/l 
servers with service rates /ii and fi2, and an arbitrary arrival process, the distribution of the 
departure process is unchanged if /ii and ^2 are exchanged. 

This interchangeability result was subsequently proved in a number of different ways, for 
example in [3], [IH], and [H]. The coupling proof given by Tsoucas and Walrand in ^\ is 
important for our purposes, since we can use their approach to extend the interchangeability 
result to multi-type queues. 

Before developing the general batch queueing model, we begin in Section [2] by giving a 
guide to the main results and methods of proof in the particular case of the continuous- 
time ./M/l queue. Since this model is already rather familiar, we give an informal account 
without introducing too much notation. (Everything is developed rigorously in later sections). 
In addition, certain aspects are simpler in the ./M/l case; for example, the service process 
has only one parameter, so all such service processes are interchangeable, and given any 
vector of arrival rates Ai, . . . , A„j (corresponding to customers of classes 1, . . . , m respectively), 
there is a unique fixed-point arrival process which is common to all ./Af(/i)/l queues for 
fi> Xi -\ h Am. 

Our general model is introduced in Section |3| which describes the set-up of a discrete- 
time batch queue. Multi-class systems are introduced in Section [4] The Bernoulli-geometric 
distribution, and corresponding queueing servers, are described in Section [5] 

Interchangeability results are given in Section [6| These extend the results for one-type 
./M/l queues described above, to cover multi-type systems and to the more general queueing 
server model. In Section [Tj we give the construction of multi-type fixed points, and prove the 
fixed-point property using the interchangeability results. The main result is given in Theorem 



7.1 (the corresponding results in the -/M/l case are Theorem |2.1| and Theorem 2.2 1. The 



proof of the interchangeability result itself is given in Section [8] 

In Section [9] we give examples of the application of the results to several of the particular 
queueing systems described above. The final example is that of the Brownian queue. Here 
the lower-priority work in the fixed-point process corresponds to the local-time process of a 
reflecting Brownian motion; this process is non-decreasing and continuous but is constant 
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except on a set of measure 0. This is an extreme case of the "clustering of lower-priority 
customers" referred to above. 



The connections with interacting particle systems are discussed in Section 10 The fixed 
points for M/M/1 servers in discrete time and in continuous time correspond to stationary 
distributions for multi-type versions of the TASEP jlSj and of Hammersley's process [H], 
respectively. Time in the queueing systems corresponds to space in the particle systems; 
with this identification, questions of fixed points for queues and stationary distributions for 
particle systems are closely analogous. 

Finally in Section 11 we mention a limit as the number of classes goes to infinity, with the 
density of each class going to 0. In this limit, the class-label of each customer becomes, for 
example, a real number in [0, 1]. This gives another illustration of the clustering phenomenon; 
although a priori, any given label has probability of occurring, nonetheless for any realisation 
of the process, each label that does occur will occur infinitely often with probability 1. 



2 Continuous-time ./M/1 case 

2.1 -/M/l queueing servers and Burke's theorem 

Let A and S be independent Poisson processes, of rate A and /i respectively, with X < fj.. We 
can use these processes to define an M/M/1 queue with arrival rate A and service rate /i. 

Arrivals occur at points of the process A; at these points the queue-length increases by 1. 
At a point of the process S, a departure occurs and the queue- length decreases by 1, unless 
the queue-length is already 0, in which case it stays the same and we say that an unused 
service has occurred. 

We write D for the process of departures and U for the process of unused services. 

Perhaps the most famous result in queueing theory is Burke 's theorem, which states that 
the departure process D is itself a Poisson process of rate A. We may regard the service 
process S as an operator (a "./M/l queue", or a "./M(^)/l queue" if we want to emphasise 
the service rate) which maps the distribution of an arrival process to the distribution of a 
service process. In this sense, Burke's theorem [7] says that a Poisson process of rate A is a 
fixed point for a ./M(/i)/l queue whenever fj, > X. 

In [4] , Anantharam showed that in fact Poisson processes are the only such ergodic fixed 
points. Mountford and Prabhakar showed further that these fixed points are attractive; 
if one starts with any ergodic process of rate A and applies the "./M(/i)/l queueing operator" 
repeatedly, the sequence of distributions obtained converges weakly to a Poisson process of 
rate A. 



2.2 Multi-class queues 

In this paper we consider fixed-points in the context of multi- class queues. 

Consider again a queue whose service process is a Poisson process of rate The queue 
may now contain several types of customer, say types 1, 2, . . . , m. Each arrival to the system is 
of one of these types. Customers of type 1 ("first-class customers") have the highest priority, 
followed by those of type 2 ( "second-class customers" ) and so on. When a service occurs in 
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the queue, if there are any customers present, then a customer with the highest priority out 
of those present departs from the system. Hence each departure from the system also has 
a type. Again one may have "unused services" , when the queue is completely empty at the 
time of an event in the service process. 

An rn-type queue may be seen as a coupling of m one-type queues, which share the same 
service process. Given r with 1 < r < m, consider only customers of classes 1 up to r, ignoring 
the differences between these customers. Since any such customer has priority any customer 
whose class is higher than r, the process obtained behaves exactly as a one-type queue. 

2.3 Multi-class fixed points 

The distribution of a multi-type arrival process is said to be a fixed point if the process of 
departures from the system has the same distribution as the process of arrivals. (We assume 
that arrivals and services are independent). 

We will consider arrival processes which are stationary and crgodic, in which case for each 
TO there is a deterministic long-run intensity of arrivals of customers of type m. 

A coupling argument analogous to that used by Mountford and Prabakhar in [23] can be 
used to show that for any Ai, . . . , Am with Ai -I- • • • -I- Xm < there is a unique stationary 
and ergodic TO-type arrival process with intensity Ar of customers of type r which is a fixed 
point for the ./M{^)/1 queue. 

As observed above, the system comprising only customers of types 1 up to r can be seen 
as a single one-type queue (for each 1 < r < to). Since the only one- type fixed points are 
Poisson processes, this shows that in any multi-type fixed point, the combined process of 
customers of types 1, . . . , r must be a Poisson process of rate Ai -|- • • • -I- A^. 

2.4 Construction of a 2-type fixed point 

First we describe how to construct such a two-type fixed point, with intensities Ai and A2. 

Consider the output process {D, U) of an M/M/1 queue as described above, with arrival 
rate Ai and service rate Ai -I- A2. We regard [D, U) as a two-type process, in which first-class 
customers occur at the points of D and second-class customers occur at the points of U . 
Write F2 = i^2,Ai,A2 foi' the distribution of this process. 

Theorem 2.1 The distribution -F'2,Ai,A2 '■s a fixed point for a ./M/l queue with service rate 
fi, for any ^ > Ai + A2. 

Note that indeed the sub-process of first-class customers is a Poisson process of rate Ai 
(by Burke's theorem) and that the combined process of first-class and second-class customers 
is a Poisson process of rate Ai -I- A2 (since D + U = S). 

However, note also that the process U of second-class customers is not a Poisson process. 
Rather, the second-class customers tend to cluster. For example, fix Ai and suppose that A2 
is very small. The process {D, U) is obtained as the output of a queue whose arrival rate is Ai 
and whose service rate is Ai -I- A2, so that it is operating very close to capacity. The long-run 
rate of unused services is only A2. However, if we observe an unused service, we know that 
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Figure 2.1: The process F2 is constructed using the departure and unused service processes 
of a one-type queue whose arrivals are Poisson rate Ai and whose services are Poisson rate 
Ai + A2. The process ^3 is then constructed using the departure and unused service processes 
of a two-type queue whose arrivals have distribution F2 and whose services are Poisson rate 
Ai + A2 + A3. 

the queue is currently empty; hence just after an unused service, the instantaneous rate at 
which another unused service occurs is Ai + A2 which may be much larger. 



2.5 Recursive construction of multi-type fixed points 



Now we show how to construct fixed points with a larger number of classes recursively. Fix 
Ai, . . . , Am- Suppose we have already constructed a distribution of an (m — l)-type process 
Fm-i, with intensity A^ of ith-class customers for 1 < i < m — 1, which is a fixed point for a 
./M(/x)/l queueing server (whenever > Ai + • • • + Am_i). Now consider a ./Af(Ai + • • • -I- 
Am)/1 queue whose arrival process has distribution F^-i- Write (-Di, D2, ■ ■ ■ , D^^-i, U) for 
the output process of this queue, comprising departures of types 1, 2, . . . , m — 1 and unused 

services. Now identify points of U as mth-type customers, and write F„i — Fm,\i a„ for 

the distribution of the m-type process obtained. 



Theorem 2.2 The distribution Fm.x^ 
rate /i, whenever /i > Ai + • • • + A,™. 



is a fixed point for a ./M /I queue with service 



The construction of F2 and F3 is illustrated in Figure 2.1 



Note that if we take the a process with distribution F,„ and ignore the mth-class cus- 
tomers, we obtain a process with distribution i^m-i- (This had to be the case, because of the 
uniqueness of the fixed-points with given intensitites, and because the ignoring the mth class 
of an TO-class system simply gives an m — 1-class system) . 



Theorem 2.2 was already proved in jj^; it emerged as a corollary of a multi-type tandem 
queue construction which was used to give the stationary distribution of a multi-type version 
of an interacting particle system called "Hammersley's process" . In this paper we give a much 
more direct proof, using properties of interchangeability of queueing servers. 
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2.6 Tandems and interchangeability 

In 1979 Richard Weber |30) proved an interchangeability property for ./M/1 queueing servers. 

Consider two independent ./M/l queueing servers in tandem, with service rates /^i and 
/i2- The first queue has some arrival process A, with an arbitrary distribution (for example, 
it could even be deterministic), which is independent of the service processes. By "tandem" 
we mean that a customer leaving the first queue immediately joins the second queue; the 
departure process from the first server is the arrival process of the second. 

The result of [30] is that the law of the departure process from the system (that is, the 
departure process from the second queue), is unchanged if /ii and /i2 are exchanged. 

Given a queue with arrival process A and service process S, we write D{A^ S) for the 
departure process from the queue. So in a system with two queues in series, with arrival 
process A and service process Si and S2 at the first and second queues respectively, the 
process of departures from the system is D{D{A,Si),S2)- Then the result of [30] can be 
written as follows: 

Theorem 2.3 Let Si and S2 be independent Poisson processes of rates fXi and ^2 respectively. 
Then for any arrival process A, 

D{D{A, Si),S2)=D{D{A,S2),Si). 

Thus, for any arrival process, the order of the queues does not affect the law of the output 
of the system. By induction, the result extends easily to tandems containing more than two 
queues. 

Alternative proofs of this interchangeability result were subsequently given by Anantharam 
[3] , by Lehtonen [19] , and by Tsoucas and Walrand [29] . The proof by Tsoucas and Walrand 
in [3S] is particularly important for our purposes since it allows an extension of the result 
which will also apply to multi-class systems. Their methods can be used to give the following 
result. 

Theorem 2.4 Let Si and S2 be independent Poisson processes of rates fii and fi2- There is 
a coupling of Si and S2 with two further processes Si and S2 such that: 

(i) (Si,S2) = (6*2,51); and 

(ii) If A is any arrival process, then D{D{A, Si), S2) = D[D{A, Si), S2) ■ 

The key point is that we can couple a system with service rates /ii , /i2 and another system 
with service rates /i2,Mi in such a way that for any arrival process, the departure processes 
from the two systems are the same. (Weber's result already implies the existence of such a 
coupling for any fixed arrival process, but not necessarily that the same coupling works for 
all arrival processes simultaneously). 

Above we explained how a queue with m types of customer can be regarded as a coupling 
of m one-type systems, with the same service process but with different arrival processes. 
(In the rth system, we only look at customers whose class in the original system is between 



1 and r, and neglect all the customers whose class is greater than r). Since Theorem 2.4 
gives a coupling which yields interchangeability simultaneously for all arrival processes, it 
immediately implies interchangeability for any multi-type arrival process: 
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Corollary 2.5 Let Si and S2 be independent Poisson processes of rates fj,i and ^2- Then for 
any multi-type arrival process A, 

D{D{A, Si),S2) = D{D{A, ^2), ^1). 



The coupling described in Theorem 2.4 is quite easy to construct. (The construction in 
|29j is essentially the same although not expressed in quite the same way). Assume without 
loss of generality that /ii < ^2- Now consider a queue whose arrival process is and whose 
service process is 82- Let D — D{Si, S2) be the departure process of this queue, and let U 
be the process of unused services (so that S2 = D + U). 

Now let S2 — D, and let Si = Si + U. Loosely speaking, we have transferred the unused 
service events from one process to the other. 

The fact that {S2,Si) has the same distribution as (6*1, is an extension of Burke's 
Theorem, and can be proved using simple reversibility arguments (see for example Theorem 
3 of 



one can 



I). To prove the intcrchangeability statement in part (ii) of Theorem 2.4 
verify pathwise that, whatever the arrival process is, transferring "unused services" in this 
way cannot affect the overall departure process. 



A full proof of this result in the more general setting is given below (Theorem 6.1 ). The 
construction in the general setting is slightly more complicated; again we remove the process 
of "unused services" from 6*2 to obtain 82, but now we do not add precisely the same process 
to 5*1 to obtain Si, but instead something more like an independent copy of the process. 
However, the same ideas of reversibility are again the key part of the proof. 



2.7 Proving multi-class fixed points using interchangeability 

Still in the context of the -/M/l queue, we finally indicate how the interchangeability result 
of Corollary 1 2. 5 1 can be used to prove that the processes defined above are indeed fixed 
points. The arguments in this section arc informal but can be made fully rigorous, as they 



are below in the proof of Theorem 7.1 We will show the argument for m — 2; the case m > 2 
follows by induction in a natural way, and can be found in the general setting below. 

Recall that F2 = i^2,Ai,A2 is the distribution of a two-type process which can be obtained 
by taking the departure and unused service processes {D,U) from an M/M/1 queue whose 
arrival rate is Ai and whose service rate is A2. So the first-class customers occur as a Poisson 
process of rate Ai, and the first-class and second-class customers together occur as a Poisson 
process of rate Ai -f A2. 

Let /i > Ai + A2, and consider feeding the arrival process F2 into a system of two queues 
in series, with service rates Ai -|- A2 and /x. The two possible orderings of the two queues can 
be visualized as follows: 



F2 
F2 







s^ 









H 
K 



For example, G is the distribution of the departure process of the first queue with service rate 
^1 + A2, and is then itself used as the arrival process for the second queue with service rate ^. 
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Note that the total arrival rate in the process F2 is Ai + A2. But the first queue in the 
first system only has service rate Ai + A2, so it is saturated; every service is used by some 
customer. Since the arrival process of first-class customers is a Poisson process of rate Ai, the 
first-class departures in G occur as the departures from an Af(Ai)/M(Ai + A2)/l queue; the 
second-class customers fill up all the other services which are unused by first-class departures. 
Hence this two-type process has the distribution of the departure and unused service process 
{D, U) from an Af(Ai)/M(Ai + Aa)/! queue, and we have G = F2. 

The same argument shows also that K = F2, since the process J has total intensity A1+A2, 
and (by Burke's theorem) its first-class customers occur as a Poisson process of rate Ai. 

Finally, the interchangeability result in |2.5| tells us that H = K, so that also H = F2. 
Since we have G = H = F2, we can concentrate only on the second queue in the first line, to 
obtain 

F2 — > Sf^ — > F2 
so that indeed F2 is a fixed point for the ./M(/i)/l queue as desired. 



3 Batch queueing model 

We now move to the more general queueing model. We begin by defining a model of a batch 
queue in discrete-time. By taking particular values, or appropriate limits, this model will 
cover various interesting cases, including discrete-time and continuous-time M/M/l queues. 

The batch queue is driven by an arrival process (A(n),neZ) and a service process 
{S{n),ne Z). 

At time-slot n G Z, A{n) customers arrive at the queue. Then service is available for S(n) 
customers; if the queue- length is at least S{n), then S{n) customers are served, while if the 
queue length is less than S{n) then all the customers are served. 

We define various processes as functions of the basic data A and S of the queue. 

Let X{n) be the queue length after the service S{n~ 1), before the arrival A{n). Formally, 
we define the process X by 

n-l 

X{n) = sup J2 (Mr) - S{r)) (3.1) 

(where a sum from 71 to n — 1 is understood to be 0) . In accordance with the description of 
the queue above, we have the basic recurrences 

X{n+l) = [X{n)+A{n)-S{n)]^ (3.2) 

(where [a;]-|- denotes max{a:;,0}). 

Let D(n) be the number departing from the queue at the time of the service S{n). So 

D{n) = mm{X{n) + A{n),S{n)) 
= X{n) + A{n) ~ X{n + 1). 

Finally let U{n) = S{n) — D{n) be the unused service at the time of the service S{n). 
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Figure 3.1: The evolution of the queue with batch services and arrivals 



See Figure 3.1 for a representation of the evolution of the queue along with its inputs and 
outputs. 

Since X, D, U are all functions of the data A and S', we sometimes write D = _D(yl, 5), 
and so on. 

Note that we allow the possibility that X(n) = oo. Indeed, we don't impose any stability 
condition on the queue; so, for example, if the average rate of arrivals exceeds the average rate 
of service, then the queue will become saturated. In fact, the following simple observation 
will be useful later: 



Lemma 3.1 Suppose the arrival and service processes are independent, and An are i.i.d. 
with mean X and Sn are i.i.d. with mean /i < A. Then with probability 1, 

Xn = OO for all n, 
Dn = Sn for all n, 
Un = for all n. 



The lemma follows immediately from the definition (3.1), since the random walk whose in- 
crement at step r is A{—r) — S{—r) is either recurrent (if n = X) oi escapes to -|-cxd with 
probability 1 (if /i < A); in either case it attains arbitrarily high values. 

If the A{n) and S{n) take integer values, then it is natural to talk in terms of "number of 
customers" arriving, or departing, or in the queue, and so on. However we will also consider 
cases where the values are more general, in which case one could talk of "amount of work" 
rather than "number of customers" . 



4 Multi-class queues 

We will now define a multi-class batch queue. The system can now contain different types 
of customers (or work) with different priorities. When service occurs at the queue, it is first 
available to first-class customers. If there is more service available than there are first-class 
customers present in the queue, the remaining service (unused by the first-class customers) is 
then offered to customers of lower class, starting with second-class customers, then third-class 
and so on. 

For example, suppose that at the start of time-slot n, there are 7 customers in the queue, 
of whom 3 are first-class, 1 is third-class and 3 are fourth-class. Suppose that 5 units of 
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service are available at time n. Then the departures at time n will be 3 first-class customers, 
1 third-class customer and 1 fourth-class customer, leaving 2 fourth-class customers remaining 
in the queue. 

Let TO be the total number of classes. We will have a collection of arrival processes 
A^^ , A^^, . . . , A^"', where A^'^{n) is the number of rth-class customers arriving at time n. 
We will also denote A-^ — X]i=i for r = 0, 1, . . . , to. 
Let S be the service process of the queue. 

Similarly we will write D^^{n) for the number of rth-class customers departing at time 
n, and X^'^[n) for the number of rth-class customers present in the queue at the beginning 
of time-slot n. Write also Z?^'^ = XiLi D=^ and X^'' = Y.\=i 

There are two natural ways to construct the multi-class queue (which are equivalent). 
One way would be to look at the queueing process of rth-class customers for each r. This 
is a queue with arrival process and service process S — D~^\ the services available to 
rth-class customers are those that have not been used by any higher-priority customer. 

Alternatively, we will consider, for each r, the queueing process of customers of classes 
1, 2, . . . , r combined. This is a queue with arrival process A-"^ and service process S. So in 
particular we define 



This second description turns out to be more useful since it describes the m-class queue 
as a coupling of to single-class queues, each with the same service process. Once we come to 
consider interchangeability of queues, the fact that we are working with a common service 
process becomes crucial. 

5 Bernoulli-geometric distribution 

We define a Bernoulli- geometric distribution, with parameters p and a. A random vari- 
able with this distribution has the distribution of the product of two independent random 
variables, one with Ber(p) distribution and the other with Geom(Q;) distribution. That is, 
A ^ Ber(p)Geom(a) if 



We have E A = p/a. 

We will consider a queue where A{n) is an i.i.d. sequence with A{n) ^ Ber(p)Geom(a), 
and S{n) is an i.i.d. sequence (independent of A[n)) with S{n) ~ Ber((7)Geom(/3). (We will 
say that A and S are Bernoulli- geometric processes). Queues of this type are investigated in 
|21j . in particular regarding their reversibility properties. 

For stability we will assume that ¥.A{n) < ES{n), i.e. that p/a < q/l3, and further we 
assume that 



D 



D(A^^s•), 



X^'' = X{A^'',S). 




c{p,a) = c{q,l3), 



(5.1) 
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where for p,a € [0, 1] we define 



cip,a)^-^- . (5.2) 

1 — p 1 — a 



When the arrival and service processes share a value of the parameter c in this way, a 
version of Burke's theorem holds. [21 . We use an asterisk to denote the time reversal of a 
process, so that A*{n) = A{—n). The following result is included in Theorem 4.1 of pT| : 

Theorem 5.1 Under the assumptions above, the departure process {D{n),n £ 1,) has the 
same law as the arrival process (A{n),n G Z). Moreover, the queue is reversible in the sense 
that 

{A{n),D{n),neZ) = {D* {n), A* {n),n e Z) . 

Note that for given q, (3 and A, there is precisely one choice of the parameters p and a 
such that p/a = X and c{p,a) = c{q,P). So for a Bernoulli-Geometric service process with 
given parameters q and /3, Theorem |5.1| gives a one-parameter family of fixed-point arrival 
processes, indexed by the arrival intensity p/a. 

In the next section we see further that two service processes are interchangeable precisely 
when they fall within the same one-parameter family (that is, when they share a value of the 
parameter c). 



6 Interchangeability 

We consider two queues in tandem. Suppose we are given an arrival process Ai and service 
processes Si, 82- From this we construct a system of two queues. The first queue has arrival 
process Ai and service process Si. The second queue has service process S2 and arrival 
process A2 defined by A2 = D{Ai,Si). (The departure process from the first queue becomes 
the arrival process from the second queue). 

Note that a departure from queue 1 at time n becomes an arrival at queue 2 at the same 
time-slot n (and may indeed depart from queue 2 at time n also). 

The following result says that two service processes, which fall within the same one- 
parameter family mentioned above, are interchangeable: 

Theorem 6.1 Suppose that the processes Si and S2 are independent of each other and each 
are i.i.d., with Si{n) ^ Ber(j))Geom{a) and S2{n) ^ Ber{q)Geom{j3). Suppose further that 



condition (5.1) holds 



Then there is a coupling 0/(51,5*2) with another pair of processes (5i,52) such that: 

(i) (5i,52) = (5*2, 5i); 

(ii) If A is any arrival process, then D2 = D2 where D2 = D(D{A, Si), S2) and D2 = 
D{D{A,Si),S2). 

Note that in particular it follows that for any A, the distributions of D{D{A, Si), S2) and 
D{D{A, S2), Si) are the same. The result gives something rather stronger; namely a coupling 
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of the two pairs of service processes such that, simultaneously for all arrival processes, the 
output of the system is the same for both pairs. 

Recall that a multitype queue can be seen as a coupling of several queues, with differ- 



ent arrival processes but the same service process. Because Theorem 6.1 gives a coupling 
which works for all arrival processes simultaneously, we immediately have as a corollary an 
interchangeability result for multitype queues. 



Theorem 6.2 Suppose the conditions of Theorem \6.1\ hold, and let A be a multitype arrival 
process. Then D{D{A, Si), S2) ^ D{D{A, S2), Si) . 



7 Multitype fixed points 

In this section we will define a series of distributions of multi-class processes, denoted by 

-^m,c, Ai , . . ., A„i ■ 

^m.c,Ai,...,A„ will be the distribution of an m-type process, with intensity of jth-class 
customers equal to Xj. We will show that it is a fixed point arrival process for a Ber(q,/3) 
service process, for all q and /3 such that q/P > Ai + • • • + A„i and c{q, /3) — c (where the 



function c is defined at (5.2)). When c and {A^} are fixed, we may abbreviate and write 
simply Fm- The distributions are defined recursively, using F^-i to construct F^. 

-Fi is the distribution of a 1-type process; simply a Ber(p)Geom(Q!) process where p and 
a are chosen so that p/a = Xi and c{p, a) = c. 

Suppose we have constructed Fm-i- Now Fm is constructed as follows. We consider a 
queue whose service process is Ber(gm)Geom(/3m), where qm and (3m are chosen so that the 
total service intensity qm/Pm is equal to Ai + • • • + A^, and so that c{qjn, I3m) = c. (As 
observed in Section [5] there is a unique such choice). 

The arrival process to the queue has distribution F„i-i (and is independent of the service 
process). This leads to an (to — l)-type departure process. We will extend the departure 
process to an TO-type process by replacing the unused service in the queue by customers of 
type m. Let Fm be the distribution of the TO-type process obtained in this way. 

Note that a consequence of the construction and of the fixed point result is that the 
restriction of a process with distribution F^ to its first to — 1 types gives a process with 
distribution Fm~i- 

Theorem 7.1 Let S be a Ber{q)Geom(j3) service process, with q/P = fj, and c{q,f3) = c. 
Then for all {Xj} with Ai + • • • + A„i < fj,, the distribution Fm,c,Ai,...,A„ is a fixed point for 
the process S. 

That is, suppose A is an m-type arrival process with distribution -F'm,c,Ai,...,A„ o,^d is 
independent of S . If ID — D{A, S) is the m-type departure process, then D = A. 

Proof: The case to = 1 is the version of Burke's theorem given in Theorem |5 . 1 1 above . The 



arrival process is a one-type Ber(p)Geom(a) process and (5.1) holds, so indeed the departure 
process has the same law as the arrival process. 

Now let to > 2 and suppose that the result holds as claimed for F„i-i- We wish to 
establish the result for Fm- We will use the interchangeability result in Theorem |6.2| 
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We consider two independent service processes and Sx-^^ |.a„, whch are Bernoulli- 
geometric processes with parameters Ber(r, 7) and Ber{qm, Pm) respectively. We choose the 
parameters such that: 

• ''/t = M sind c(r, 7) — c; 

• q-m/Pni = Ai H h Am and c{q,n,Pm) = c. 

We will consider an m-type arrival process A distributed according to Fm- We con- 
sider the effect of feeding this arrival process into two tandem systems, comprising the ser- 



vice processes S'^ and S\^^ |-a„ in their two possible orders. Theorem 6.2 tells us that 

£'(-D(A, S"^), ^Ai+.-.+A,,,) and D{D{A,Sx^^ hA„J,'S'^) have the same distribution. To make 

use of this result we first need two simple properties. The first concerns projections of multi- 
type processes onto processes with fewer types, and the second concerns the behaviour of a 
queue which is saturated (i.e. whose service rate is only just sufficient to serve all the arriving 
customers). 

Claim 1: Suppose A has distribution i^„i, and let G,„ be the distribution of D{A^S^). 
Then the restrictions of F.„i and Gm to the first m—1 coordinates have the same distribution, 
and the combined process of all m types under Gm has the same distribution as S'a^h hA„, • 

Proof of Claim 1: Restricting Fm to its first m—1 coordinates gives the distribution 
Fm-i- Since Fm-i is an (m — l)-type fixed point for the service process 5^, the first m—1 
types under Gm do indeed have distribution Fm-i also. 

Finally, from the defintion of Fm, the combined process of all m types under Fm is a 
Ber(gm)Geom(/3m) process. The service process is a Ber(r)Geom(7) process. By assump- 



tion c{qm,l3m) — c{r,j) — c. So condition (5.1) is satisfied for the queue, and we can apply 



the 1-type fixed point result in Theorem |5.1| to show that the distribution of the combined 
process of all m types in the departure process has the same distribution as the combined 

process of all m types in the arrival process, namely Sx-^-\ • This completes the argument 

for Claim 1. 

Claim 2: Suppose A, restricted to its first m—1 coordinates, has distribution Fm-i- 
Suppose also that the combined process of all m types of arrival in A is an i.i.d. process with 
intensity Ai + ■ ■ • + Am. Then D{A^ Sxi-\ i-a„, ) has distribution Fm- 

Proof of Claim 2: Recall that Fm is the distribution obtained by passing an (m — l)-type 

arrival process with distribution Fm-i through a queue with service process Sxi-\ i-a„, i and 

putting customers of type m in place of all unused service. 

The first m — 1 components of A indeed have distribution Fm-i- Certainly the mth 
component of the departure process is a subset of the unused service from the m — 1 first 
types. To show that the departure process has distribution Fm, it remains to show that all 
the service unused by the first m—1 types is used by customers of type m; that is, that if we 
look at all m types combined, then there is no unused service. 

The arrival process of all customers combined is an i.i.d. process with rate Ai + • • • + Xm] 
the same is true of the service process, and the arrival and service processes are independent. 
Hence by Lemma [XT] there is no unused service in the queue as desired. So indeed the output 
process has distribution Fm, and we have established Claim 2 as required. 
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Now we put together Claims 1 and 2 to complete the proof. We consider the two possible 

orderings of the service processes S'^ and S\-^^ The two systems that arise can be 

illustrated as follows: 



F„ 

Frr 



Su 



G„ 
Jill 



SxiA hA„ 



H„ 



By Claim 1, Cm satisfies the conclusion of Claim 1, which is also the assumption on the 
arrival process of Claim 2. Then by Claim 2, Hm = Fm- 

Certainly F,„ itself also satisfies the assumption of Claim 2, so by Claim 2, J™ = Fm- 
Theorem |6. 1| says that = K^- 

Then the second box of the second line down shows us that if we pass J„i through we 
get Km- 



But J„i = Km = Fm- So indeed Fm is a fixed point for as desired. 



□ 



8 Proof of interchangeability result 



In this section we will prove Theorem |6.1| We first need two simple lemmas about queues in 
tandem. 

Lemma 8.1 Consider a system of two queues in tandem with arrival process A and service 
processes Si and 5*2, all taking integer values- For s G Z. define oo) to he the arrival 

process truncated before time —s. That is, 



if n < —s, 

A{n) if n > —s. 



Define D2 = D{D{A, Si), S2) and Dj^-* — D{D{A\(^_g^^j, Si), S2). For any n, if s is large 
enough then D2{n) = D\^\n). 

Proof: We compare the original system and the system with the arrival process before time 
—s. Write Di and D2 for the departure processes from the first and second queues in the 
original system, and similarly D\ , D)^ for the processes in the system with the truncated 
arrival process. 

Take any n From the definitions in Section |3j we have 

Di{n) = min [Si(n),A{n) + sup ^ \A{r) - Si{r)\ \ , 

\ "<« r=u ) 

D['\n) = min (si{n),A{n) + sup ^ [A{r) ~ 5i(r)] ) . 



-s<u<n 
— — r—u 



Certainly D[^\n) is increasing in s, and D[^\n) < Di{n) for all s. 

We first observe that D[^\n) — Di{n) for all large enough s. Suppose the sup in the first 
line is finite. Then (since the variables all take integer values) it is attained for some u = u*, 
and we have Di{n) — D[^\n) for all s > —u*. If instead the sup in the first line is infinite, 
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then by taking s large we may make the sup in the second hne as large as desired; if s is large 
enough that the second term in the min exceeds 5*1(71), then we have Di(n) — -D[^''(n), as 
required. 

Now we have the corresponding expressions for the departure process from the second 
queue (whose arrival process is Di and whose service process is 5*2): 

D2{n) = min [s2{n),D^{n) + sup^ [7^i(r) - S2{r)] \ , 

V r=u J 

D['\n) = min ( S2{n), D['\n) + sup V \D^^\r) - 52 (r)!) . 

Since D''-^\n) is increasing in s and bounded above by Diin), we also have that D2^\n) is 
increasing in s and bounded above by D2{n). 

Suppose that the sup in the first line is finite. Then as before it is attained at some u* . 
Now if we take s large enough that D''^\r) — Di{r) for all r with u* < r < n, then we see 
from the second line that indeed D'"2\n) > D2{n), and hence in fact D2^\n) — D2{n) as 
required. A similar argument applies if the sup in the first line is infinite; by taking s large 
enough, we can make the sup in the second line as large as desired. □ 

Lemma 8.2 Consider a tandem of two queues, with arrival process A and service processes 
Si and 82- Suppose that A{n) ~ /or all 71 < 0. Let D2 be the departure process from the 
second queue. Then for all t>\, 

t-l (ux-l U2-I t-1 ^ 

E^(^) = o<„-4<, E + E ^i(^) + E S.ir) . (8.1) 

r— k r—0 r—ui r—U2 ) 

Proof: An equivalent result for a related queueing model was given in |28j , and similar 



properties appear in many places in the literature. Formula 8.1 is essentially a special case 
of equation (10) of [27]. 

Let Dx and D2 be the departure processes from the first and second queues, and Xi and 
X2 be the queue-length processes at the first and second queues. 

Since An — Q for all negative n, the queues start empty at time 0: A'i(O) — ^"2(0) = 0. So 
the total number of departures from the system (that is, from the second queue) before time 
t is given by the total number of arrivals during that time, minus the queue-lengths at time 
t: 

t-1 t-1 

^ D2{r) = ^(^) - - ^2(t). (8.2) 

r=0 r=0 

Similarly, considering the evolution of the first queue alone between times u and t, we 
have that for any u < t, 

t-1 t-1 

^i?i(r) = ^A(r)-Xi(i) + Xi(«). 

r—u r—u 

Now X2 is the queue-length of a queue with arrivals Di and services 52. So from the 
definition (3.1), we have 



Xi{t) + X2(t) = Xi{t) + sup J2 [^i(^) - '^2(r 



^— r—U2 
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+ sup i V A{r) - Xi{t) + X^{u) - V 52 (r) 



^r—U2 r—U2 
t-1 "1 



snp\x,{u)+y] [A{r)-S2{r) 

U2<t I 



r=U2 



{U2-1 t-1 
sup £ [A(r) - 5i(r)] + J2 [Mr) - S2{r)] 
i_ z J-— r—U2 ^ 

{t-l t-l \ 

Y,A{r)-Y,S,{r)-Y,S2{r)\. 
r—ui r—ui r—U2 ) 



Since A„ = for negative n, the sup will be attained for some u\_ > 0. Then using (8.2 ) gives 
(8.1 1 as desired. □ 
Now we proceed to prove Theorem 6.1 To do this, we need to construct Si, 82,81 and 
82 on the same probability space, in such a way that: 

(i) {81,82) = (82,81) 

(ii) For all A, D{D{A, 81), 82) - D{D{A, 8^), 82). 

We assume that ^2 has higher intensity than 81. (If not, then swap them around. If the 
intensity is the same, then the assumptions ensure that they have the same distribution.) 

The construction is as follows. We think of a queue with arrival process and service 
process 5*2. We can decompose 82 into -0(51, 5*2) and U{8i,82) - that is, 5*2 = 0(81,82) + 
U{8u82). 

We let ^2 — D(8^ , 82) - So we have obtained 82 from 82 by removing the unused service 



process. By Theorem 5.1 we have that (5*1, 82) = {82, 8^). (Recall that the asterisk denotes 



the reverse of a process). Knowing this, we can extend the probability space to include a 
random variable Si, with its time-reversal Si, such that {81,82, 82) = {S2, 81,81). 

In particular this implies that {81,82) — {82,81). But the distribution of {81,82) is 
invariant under time reversal, so this definition gives us (i) above as desired. 

(Note that since {81,82,82) = {8^ , 8*1 , S*i) , we have in particular that 8*1 = D{S^,S*i), 
or equivalently that 8^ + U{S2, 8^) — 8^. So we can think of as being obtained from 
by adding an unused service process, for a queue operating in reverse time). 



Now we need to verify (ii) above. Using Lemma 8.1 it will be sufficient to verify it for 



every arrival process which is up to some finite time. Without loss of generality, suppose 



that A{n) = for aU n < 0. Let D = D{D{A, 81), 82)). From Lemma |0| we have the 
following representation of the total number of departures from the tandem of two queues 
from time up to time t — 1: 

t-1 (ui-l U2-I t-1 ~j 

E^w - o<„-4<, E + E + E ^^(^) ■ 

r— k r—0 r—ui r—U2 ) 

It will be enough to show that for any A, this quantity is unchanged if we replace (^i, 82) by 
(5*1, 82). For this, in turn it's enough that for all s and t, the quantity 
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is unchanged if we replace (5*1, by (5*1, 6*2). 

We consider the queue with arrival process Si and service process 5*2 , along with its queue- 
length, unused service and departure processes X, U and D which are functions of Si and 
S'2. Note that U{n) > precisely if S'2(n) < 52(n). 

First we show that the quantity we are interested in is unchanged if we replace (>S'i,S'2) 
by (5'i, 5*2); that is: 

i-f<, (E + g S,ir)\ ^ Tf Siir) + g ^2(0} ■ (8.3) 

r—s r—u ) \ r—s r—u ) 

If there is no unused service between times and i — 1, then S'2 and ^2 agree on the whole 
interval and the equality is obvious. 

Otherwise, let n be the latest time in {s, ...,< — 1} for which there is unused service, i.e. 
for which S^in) > S2in), or equivalently that U{n) > where U = [/(Si, S'2). In particular 
there are no more arrivals than departures in the interval {s, . . . , 71} (otherwise there would 
still be customers in the queue after time n, and hence there could not have been unused 
service at time n). That is, X]r=s '^i(*) — ^"=8 ^^ii)- This implies that we can take u > n 



in minimising the RHS of (8.3). Then consider using the same u on the LHS as on the RHS. 
Since u > n and n was the last moment of unused service, the processes S2 and S2 agree after 
time n. So in fact the LHS is at least as small as the RHS. But since S2 > S2 it is also clear 
that the LHS is no smaller than the RHS. So in fact the two are the same. 
Now we want to show that we can further replace Si by Si; that is, 

>J<, (E + g S2(.)| ^ ^ inf^^ (g Si(.) + g S2(0| . (8.4) 



This is in fact precisely the same relation as (8.3 1, but now for the queue with arrival process 



S2 and service process Si- So an exactly analogous argument can be applied. 



This completes the proof of (ii) and hence of Theorem 6.1 



9 Summary and examples 

In this section we point out various interesting cases that can be obtained by taking particular 
values or appropriate limits in the general model considered up to now. 

Our general model is of a service process in discrete time, whose batches are i.i.d. with 
Ber(g)Geom(/3) distribution. Such a process has service intensity q/(3. This model has two 
parameters. For each choice {q, /3) there is a one-dimensional family of service process which 
are interchangeable with it; namely an (r, 7)-process is interchangeable with a (g, /3)-process 
whenever c(r, 7) = c(g,/3). This relation partitions the parameter space into equivalence 
classes (and each class contains precisely one set of parameters for each service intensity). 

Suppose we are given service parameters {q,P) with service intensity n = q/(3, and also 
Ai, . . . , Am such that Ai + • • • + A„i < /i. Then there is precisely one m-type fixed point arrival 
process with intensities A^ of rth-class customers for 1 < r < m. For example, the two-type 
fixed point with intensities Ai and A2 is obtained by taking the departure and unused service 
processes from a queue with arrival process A and service process S where A and S are chosen 
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so that A has intensity Ai, S has intensity Ai + A2, and both A and S are interchangeable 
with the {q, /?)-process. (As explained in the last paragraph, this uniquely determines A and 
S). 

This framework extends to various particular cases which we describe below. In some 
cases, the model has only one parameter, and any pair of processes in the class considered 
are interchangeable. 

9.1 -/M/l queue in discrete time 

If we take l3 — 1, the service process is a "Bernoulli process". The batches are i.i.d., each 
equal to with probability I ~ q and 1 with probability q, so that fi = q. Hence we have a 
one-parameter model. 

This is the discrete-time equivalent of the -/M/l server described in Section [2] A version 
of Burke's theorem for such systems was proved by Hsu and Burke in |17) : any Bernoulli 
process with lower intensity p < q is a one-type fixed point arrival process. 

Any two such Bernoulli processes are interchangeable. For any m and Ai + • • • -I- A„i < 1, 
there is an arrival process -Fm,Ai,....A,„ with intensities A; and which is a fixed for any such 
./M(/i)/l server whenever /i > Ai -(-••• -I- \m- The multi-type fixed points were constructed 
in |13j . and they can also be seen as stationary distributions for multi-type versions of the 
TASEP (totally asymmetric simple exclusion process). 

9.2 Geometric and exponential batches in discrete time 

If we take q — 1, then the service batches are i.i.d. geometric distributions. (A similar case 
arises for q — 1 — (3, except that the geometric distribution obtained starts from rather 
than from 1). Again we have a one-parameter model, and any two such geometric processes 
are interchangeable. A related result (for a particular form of the arrival process) was given 
by Draief, O'Connell and Mairesse in [5]. A version of Burke's theorem had been given by 
Bedekar and Azizoglu in [5]. 

The multi-type fixed points in this case are related to stationary distributions for a multi- 
type version of the totally asymmetric zero-range process |22| . 

Now take the limit (5 — 0, putting j3 = 5/u and rescaling work by a factor 5. In this way 
we can obtain the case where the batches are i.i.d. exponential with mean u. Again, any two 
such exponential processes are interchangeable. 

9.3 Batch queues in continuous time 

Suppose we let e — > 0, take q — ve and rescale time by a factor e. 

Then we have batches which are geometric with parameter /?, which occur at times of a 
Poisson process of rate v. 

Now we have a two-parameter distribution, and for each pair of parameters {v^ j3) there is 
a one-parameter family of processes, one for each service intensity, which are interchangeable 
with it. Namely, processes of this kind with parameters (t^, /?) and (p, 7) are interchangeable 
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1-/3 1-7- 

If we further take 5 — )■ 0, putting 15 = 5/u and rescaling worli by a factor 5, then we obtain 
the case where the batches are exponential with mean u (again occurring at the times of a 
Poisson process of rate v). Two such "Poisson-exponential" processes with parameters {v^u) 
and (/o, w) are interchangeable if v/u^p/v. 

If on the other hand we take /3 — )■ 1, the batches all have size 1 and we obtain the familiar 
./M/1 queueing server which we described in detail in Section [2] 

9.4 Brownian queues 

Queues in continuous time (as in the previous example) can be represented using a notation 
analogous to that introduced for discrete-time queues in Section |3] 

Let S't, t € M be a service process and A(, t G M be an arrival process. We interpret St — Sg 
as the amount of service offered in the interval (s, t], and similarly At — Ag as the amount of 
work arriving in (s,t]. 

Adding a constant to either of the processes A and S makes no difference, so we may take 
for example = 5*0 = 0. 

We define processes Qt, Dt and Ut: as functions of the processes At and St, by 

Qt = ^n^\{At - Ag) - {St - Sg)], 

s<t ^ J 

Dt = At + Qo-Qt, 
Ut^St-Df 

Now Dt — Da is the amount of work departing in (s,i], Ut — Us is the amount of unused 
service in (s,t], and Qt is the queue-length at time t. 

Note that under these definitions we have Uq = Dq — (again, this normalization is not 
important). 

A frequently studied example is that of the Brownian queue. Let At be a two-sided 
Brownian motion with drift A and variance 1, and let St be a two-sided Brownian motion 
with drift jj, and variance 1, where A, /x G K and X < fi. 

This model arises naturally as the scaling limit of queues in the so-called "heavy traSic" 
regime [T^, |3I1. For example, for n G N consider a discrete time queue whose arrival process 
A{r),r e Z is a Bernoulli process with rate l/2 + \/y/n and whose service process S{r),r e Z 
is a Bernoulli process with rate 1/2 + /i/^/n. We recentre and rescale by defining, for t e M, 

* Vn/2 ' * 

Then as n — )• 00, the distribution of the processes and 5^"^ converges to that of At and 
St above. 

A version of Burke's theorem holds for such a Brownian queue: namely that the processes 
Dt and At have the same distribution. This was shown by Harrison and Williams in |16j (see 
also |26j in particular, where the result is generalized in a number of ways). The process Qt 
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has the law of a stationary reflected Brownian motion with drift — (/i — A), and its stationary 
distribution is exponential with rate fi — X. The process Ut is the local-time process of Qt at 
0. Note that U is non-decreasing, and grows at average rate fj. — X, but is constant except on 
a set of measure (namely, the set of times when (5 = 0). 

Thus we may say that a Brownian queueing server (with drift fi) has Brownian motions of 
drift A < /i as fixed point arrival processes. Again this can be extended to multi-type arrival 
processes. 

Given a service process St and a two-type arrival process A^^"^ , Ap'* , the two-type departure 
process d[^\ D^^^^ is defined in a natural way by 



Suppose we wish to construct a two-type fixed point for a Brownian service process of 
drift /i, which has drifts Ai and A2 where Ai < Ai -|- A2 < ijl. To do this, consider a Brownian 
queue whose arrival process is a Brownian motion with drift Ai, and whose service process is 
a Brownian motion with drift A2. The process {D, U) of departures and unused service from 
this queue is a two-type fixed point of the required type. 

In this example we can see an extreme form of the "clustering of lower-class customers" 
already observed in the discrete case. Here the second-class work, represented by the local- 
time process U , is concentrated on a set of times of measure 0. Note also the difference in 
nature between the process D of first-class work (which is non-monotonic, and in fact has 
unbounded variation) and the process U of second-class work (which is non-decreasing) . 

Fixed points with larger numbers of classes can be constructed recursively in the same 
way as in the discrete case. The subprocesses of mth-class work in these fixed points are 
non-decreasing and singular, for each m > 2. 

One can also obtain an interchangeability result; two Brownian service processes, with 
different drifts but with the same variance, are interchangeable. 

Arrival processes obtained as the local time of a reflecting diffusion have been considered 
before in various contexts, for example in [20] and [TH]. These models were also motivated in 
part by the modelling of a second-class departure process as a local-time process. However, 
the models analysed are somewhat different; the service processes are deterministic rather 
than Brownian {St = ct) and the local-time arrival process has access to all the service 
capacity (rather than sharing the queue with a higher-priority stream). 

10 Interacting particle systems 

There is a close analogy between fixed-point arrival processes for queues and stationary dis- 
tributions for interacting particle systems. 

Consider for example the TASEP (totally asymmetric simple exclusion process). In the 
one-type version of the process, each site of Z contains either a particle or a hole. The state- 
space of the process can be written as {1, 00}^, where 1 denotes a particle and 00 denotes a 
hole. The dynamics of the process are as follows: each particle tries to jump to its left as a 
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Poisson process of rate 1, and a jump succeeds if the site to the left of the particle is empty 
(in which case the values 1 and oo are exchanged between the two sites). 

One can also consider a multi-type TASEP. The process with m types has state space 
({1, 2 . . . , m} U {oo})^. As before, each particle tries to jump to its left as a Poisson process 
of rate 1, and a jump by a particle of type r succeeds if the site to its left is occupied by a 
particle with lower priority (higher-numbered class) or is empty (in which case the values at 
the two sites are interchanged as before). As in the case of multi-type queues, this process 
can be seen in a natural way as a coupling of m single-type TASEPs. 

Stationary distributions for TASEPs with two types of particle were constructed by Der- 
rida, Janowsky, Lebowitz and Speer [5], and related construction were also given in [TT] . 
[5] and [TU]. In [T3], this construction was extended to multi-type processes, and made 
particularly explicit using a construction using queues in tandem which corresponds to the 
construction of multi-type fixed points for queues above. 

Given Ai, . . . , A™ with Ai + • • • + Am < 1, consider the fixed-point arrival process distri- 



bution -F'm,Ai,...,A„ for the discrete-time ./M/1 queue as described in Section 9.1 Now regard 
this distribution as a distribution over configurations of the m-type TASEP; an arrival of type 
r at time n corresponds to a particle of type r at site n, and an empty arrival slot at time n 
corresponds to a hole at site n. 

The result of [T3] is that -F'm,Ai,...,A„ is a stationary distribution for the m-type TASEP 
(indeed, it is the unique ergodic stationary distribution which has intensities Ar). 

Note that time in the context of the queue corresponds to space in the context of the 
particle system. 

The analogy between queues and particle systems can be made even clearer by considering 
particle systems evolving in discrete time rather than continuous time. The same distributions 
Fm are also stationary for certain versions of the multi-type TASEP in discrete time [53]. 
(Here many particles may try to jump at the same time-step, and there are various possible 
natural conventions for the order in which the jump attempts are to be processed). Now a 
queueing server, which uses a service process to transform an arrival process into a departure 
process, can be seen as analogous to a set of jump attempts at particular sites at a given time- 
step of the particle system, which transforms the particle configuration before that time-step 
into the new particle configuration after that time-step. Successive updates in the particle 
system correspond to consecutive queueing servers in a system of queues in tandem. 

The fixed-point property for the distributions F^^ can be derived from the result for the 
particle system (as was done in [13] )■ To see this, note that if we take an m-class TASEP 
and treat particles of class m as holes, then the resulting restriction of the process to the 
first TO — 1 types is an (m — l)-class TASEP. Hence restricting the stationary distribution 
Fm of the TO-class system to its first to — 1 types must give a stationary distribution of the 
(to — l)-class system, namely Fm-i- That is, Fm~i can be obtained from Fm by restricting 
to the first to — 1 types. But also, the first to — 1 types of Fm are defined as the departure 
process from a queue whose arrival process is Fm-i (the remaining TOth class is the unused 
service process). So indeed F^-i is also the distribution of the departure process from the 
queue; this is the fixed-point result we obtained earlier. 

However, such a proof via the particle system is rather indirect and perhaps less natural 
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than the approach using interchangeability above. 

In fact, in some cases, results previously obtained for stationary distributions of particle 
systems can themselves be derived from the fixed-point results for queues. For example, 
the stationary distributions for "Hammersley's process" P correspond to the fixed points 
for the ./M.l queue in continuous time, described in Section [2] [Ml [12]. The dynamics of 
Hammersley's process can be obtained by taking appropriate limits starting from a M/M/1 
queueing server in discrete time (first one lets the rates of arrival and service approach 1 
together in an appropriate way; after rescaling space, the particles in a discrete-time version 
Hammersley's process correspond to the gaps in the customer processes in the queue. Then 
the more familiar continuous-time version of Hammersley's proces can be obtained by letting 
the rate of jumps in the discrete version go to and rescaling time). In this way the result 
on stationary distributions for Hammersley's process can be seen as a consequence of the 
fixed-point results given here. (However, we don't know how to use similar methods to derive 
the results for the TASEP, for example). 

11 Continuous class-labels 

By considering the limit as the number of classes goes to infinity, and the intensity of each 
class goes to zero, one can arrive at models where the classes are indexed by real numbers. 

This could be done in the general case above, but for simplicity we discuss it in the case of 
a -/M/l server (where the service process is a Bernoulli process). We will consider an arrival 
process in which at most one customer arrives at each time slot, and a customer arriving 
at time n has a real-valued label. As before, lower-labelled customers have priority over 
higher-labelled customers. 

By taking appropriate limits in the results above, one can obtain the following result: 

Proposition 11.1 There exists a unique distribution over sequences {L{n),ne Z) e [0,1]^ 
with the following properties: 

(i) The distribution is stationary and ergodic. 

(it) i(0) - C/[0, 1]. 

(Hi) Let < A < 1. Define an arrival process as follows: if L(n) < A, then a customer with 
label L(n) arrives at time n, while if L(n) > X, then no customer arrives at time n. 
This arrival process is a fixed point for the ./M{p)/l queue for all jj, > X. 

Property (ii) is included just as a normalization (since the mechanics of the queue are un- 
changed if the labels of all the customers are transformed by an increasing function). Then 
in property (iii), the condition fi > X ensures that the queue is not saturated. 

The distribution of {L{n),n e Z) is also a stationary distribution for the TASEP. See [5] 
for an investigation of this process, including in particular an interpretation as the "speed 
process" for a multi-type TASEP started out of equilibrium. 

One interesting property of the process is a manifestation of the "clustering" effect de- 
scribed above. Although L{n) has a continuous distribution for each n, nonetheless one has 
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that for any n, P(i(0) = L{n)) > 0. (For example, P(i(0) = i(l)) = 1/6). In fact, with 
probability one there exist infinitely many n such that L(0) = L{n). Hence clustering occurs 
in the following sense: although any class-label has probability of being seen a priori, if one 
sees the label at any particular time the same label has high probability of being seen nearby, 
and will be seen infinitely often in the process. 
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